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The extended application of graphene-based electronic devices requires a bandgap opening in order to 
realize the targeted device functionality. Since the bandgap tuning of pristine graphene is limited to 
360 meV, the chemical modification of graphene is considered essential to achieve a large bandgap opening 
at the expense of electrical properties degradation. Reduced graphene oxide (RGO) has attracted significant 
interest for fabricating graphene-based semiconductors since it has several advantages over other forms of 
chemically modified graphene; such as tunable bandgap opening, decent electrical properties, and easy 
synthesis. Because of the reduced bonding nature of RGO, the role of metastable oxygen in the RGO matrix 
is recently highlighted and it may offer emerging ionic devices. In this study, we show that multi-resistivity 
RGO/n-Si diodes can be obtained by controlling the RGO thickness at a nanometer scale. This is made 
possible by (1) a metastable lattice-oxygen drift within bulk RGO and (2) electrochemical ambient hydroxyl 
(OH) formation at the RGO surface. The effect demonstrated in a p-RGO/n-Si heteroj unction diode is 
equivalent to electrochemically driven reversible electronic manipulation and therefore provides an 
important basis for the application of O bistability in RGO for chemical sensors and electrocatalysis. 



Extensive studies have been performed to pursue the application of graphene to electronic devices because of 
its superior electronic properties. Upon chemical vapor deposition on a metal 1 (e.g., Cu) or epitaxial growth 
on a SiC substrate 9 , single layer graphene (SLG) shows the highest mobility of any material with a two 
dimensional (2D) monolayer of sp 2 -bonded carbon atoms. Because of this reason, SLG is physically classified as a 
zero-bandgap Dirac material 2 . Although the electrical properties of SLG are ideal for electrodes, the utmost 
interest in graphene is to transform it to a semiconductor for achieving the various functionalities required from 
electronic device components. A representative example related to this issue is the graphene -channel transistor. 
In order to meet the modern low-power transistor specifications, the on/off current ratio should be at least > 10 6 3 . 
However, the best result on graphene transistors revealed a value of <10 2 for the on/off current ratio 4 . This 
problem with graphene transistors is essentially associated to the lack of a sufficient bandgap for device function- 
ality 5 . There have been various approaches to achieve a bandgap opening in graphene, such as bilayer graphene 
(A-B stacking) 6 , structural modifications 7 , and chemical functionalization using adsorbed molecules 8 ' 9 . Among 
these approaches, RGO has been widely researched upon since it is easy to fabricate and offers large bandgaps (e.g. 
~3 eV) 1011 . It is obvious that the in-plane electrical conductivity of graphene oxide (GO) is significantly less than 
that of pristine graphene 12 since O 2p orbitals form sp 3 bonds as a result of the breakup of sp 2 carbon bonds 11 . 
Nevertheless, RGO with a large bandgap and decent carrier mobility, as compared to GO, has been proven to be 
suitable as a 2D semiconductor for functional electronic devices such as photoluminescencent devices, and highly 
sensitive optical 12 and chemical sensors 13 . In recent researches on RGO properties, the interaction between 
graphene and oxygen has attracted significant interest. It has been well established that O -bonded graphene 
shows p-type electronic characteristics. This was easy to establish since O or water adsorption on the grapheme 
surface from air exposure for several days leads to p-type doping of graphene 14 . More interestingly, this p-type 
character is easily modified to n-type upon vacuum annealing involving O-related molecule desorption. Y. Sato 
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et al reported that O adsorption/desorption on a bilayer graphene 
surface can be controlled by a vertical gate bias in a field effect 
transistor structure; this observation has been declared as electrically 
modified charge doping effect 15 . On the basis of these precedent 
studies, it is clear that the graphene O bond has some degree of 
bistability, which can be controlled by a relatively low external energy 
source such as a low electric field and heat. Therefore, the reliable 
control of graphene O can be a strategy for realizing graphene-based 
electronic devices. It was recently proved that the use of O bistability 
in a graphene/Cu junction combined with electric-field- driven 
out-of-plane O migration provides a resistive switching operation 
across a nanometer- scale junction 16 . It is thought that O bistability 
in graphene originates from the unique chemical bonding structure 
of 2D carbon layers. In-plane electrochemical reduction and oxida- 
tion of GO has been demonstrated in a previous study 17 . A great 
advantage of O bistability control through out-of-plane O migration 
upon the application of an electric field is that the local change of O in 
the range of several atomic distances at a very low potential is not a 
"local chemical change" but equivalent to "global charge doping," 
affecting the entire electronic properties of single- or multilayered 
graphene. 

In view of the perspectives of O -bistability in graphene, as claimed 
in previous reports described above, the present study demonstrates 
that multiresistivity RGO/n-Si diodes can be obtained by controlling 
the RGO thickness at the nanometer scale. This is made possible 
by (1) the electrochemical metastable lattice-O drift within bulk 
RGO and (2) ambient hydroxyl group (OH) adsorption on the 
RGO surface. For a typical Schottky-diode structure consisting of 
p-type RGO/n-Si, a triple bias-dependent modulation of the 
forward current conductivity is observed. On the basis of electrical 
analysis under controlled air, vacuum, and nitrogen atmospheric 
conditions, it is found that the internal lattice-O ion drift and elec- 
trochemical OH-group adsorption/desorption in the RGO layer, thin 
enough for an electrically driven ion drift, act as (1) charge injection 
barrier and (2) strong p-type doping enhancer, respectively. 
Therefore, this finding provides an important basis for the applica- 
tion of O bistability toward realizing RGO -based diodes and chem- 
ical sensors. 

Results and Discussion 

Thickness -controlled RGO films were fabricated on a preheated Si 
substrate by spin coating. From the height analysis of the step edges 
of RGO by quantitative atomic force microscope, the resultant RGO 
film thicknesses were measured to be about 390 and 5.7 nm for the 
thick and thin films, respectively (Figs. la-b). Raman scattering can 
be employed to directly probe the formation of RGO. Raman G, D, 
and 2D bands of the thin RGO sample are shown in Fig. lc. Both the 
G and 2D modes shift to significantly higher frequency upon oxida- 
tion, which is a strong evidence of RGO formation. The in-depth 
chemical properties of thick and thin RGO and full GO (FGO) were 
investigated using X-ray photoelectron spectroscopy (XPS). 
Figure 2a shows C 1 s spectra with Gaussian-fitting deconvolution 
for each sample. In the C 1 s XPS spectra, C-C (-284.5 eV), C-O/C 
= O (-286.8 eV), and C-OOH (-288.7 eV) bonds are resolved by 
peak deconvolution. The peak assignment was based on a previous 
analysis of FGO films 18 . Depending on the sample, there are distinct- 
ive bonding energy shifts; compared to thick RGO, (1) thin RGO 
shows a higher bonding energy of the C-C bond, possibly originating 
from a contribution of partial C-OH bonds, and (2) FGO shows a 
higher bonding energy of the C-O bond, possibly originating from a 
contribution of partial C = O bonds with a higher O-coordination of 
the C species. In order to clarify the trend of bonding states for each 
sample, the O/C stoichiometric ratio and C 1 s bond fractions are 
plotted in Figs. 2c and 2d, respectively. As expected, the C-O/C = O 
fraction in FGO is the highest since more O is incorporated. An 
increase in the O/C stoichiometric ratio from 0.4 in thick RGO to 



0. 9 in FGO is observed. The C-OOH bond fraction is, however, 
highest at 18% for thin RGO compared to 5-6% for thick RGO 
and FGO. Because the same deposition method was applied for thin 
and thick RGO films, except that the spin coating conditions were 
varied for thickness control, the higher level of OH incorporation in 
thin rather than in thick RGO films is believed to originate from 
different surface conditions. The process conditions for thin RGO 
provide a more favorable surface to water adsorption and its sub- 
sequent dissociation into OH groups bound to C in graphene. In 
addition, the relative surface portion in the whole film, where water 
adsorption/dissociation from air or solution occurs, significantly 
differs between two samples. Typically, the surface oxygen and water 
adsorption is limited at the top several nm- thick region and this 
inherently leads to OH/O-rich surface dominant thin RGO layer 
compared to the thick RGO film where the chemistry of entire film 
is not highly affected by the surface composition. The degree of water 
adsorption may also be related to the surface roughness as a key 
factor to provide the reactive sites (i.e., epoxy group configuration) 
for O adsorption 11 . Therefore, this proves that the surface of thin 
RGO is intrinsically or thermodynamically subject to the ther- 
modynamic water dissociation and OH incorporation in RGO. In 
addition, the important implication of OH-rich RGO is strong p-type 
electronic conductivity, which affects the RGO diode structures 
significantly. 

The valence band edge XPS spectra (Fig. lb) indicate that RGO has 
a significant metallic contribution of an sp 2 C-C bonding configura- 
tion, resulting in an overlapping of the Fermi energy level (E F level, 

1. e., 0 eV binding energy) with the valence band onset, while E F is 
— 2 eV above the valence band onset in FGO. This result implies that 
an incomplete O -coverage of the surface results in partial metallic 
behavior; this finding is in agreement with an earlier report on the 
significant bandgap modulation of RGO 19 . GO can be either a p-type 
or an n-type semiconductor, depending on the environment and the 
substrate 20 ; however, p-type behavior is more common under ambi- 
ent conditions, suggesting that 0 2 /water adsorption plays a role in 
the p-type character of RGO. The considerably less mobility than 
pristine graphene as a result of sp 3 C-O bond formation is unavoid- 
able in RGO; nevertheless, this mobility drop can be mitigated by a 
large bandgap opening. 

A log I-V plot for thin and thick RGO/n-Si diodes is shown in 
Fig. 3a. The diodes show clear rectifying junction characteristics with 
rectification factors of 64 and 4 for thin and thick RGOs, respectively, 
at a bias of ± 1 V. The inset of Fig. 3a shows the structure of the RGO/ 
n-Si diode used in this study. The Ag electrode offers a good ohmic 
contact to RGO, which is in contrast to the slightly rectifying junction 
for FGO (see Supplementary Information, Section 1). Therefore, the 
rectification occurring at the RGO/n-Si diode is truly related to 
the interfacial band alignment of the RGO/n-Si junction. Although 
the forward current in the thick RGO diode seems to exhibit a normal 
behavior when compared to regular Schottky diodes, the thin RGO 
diode shows abrupt current saturation at a relatively low voltage of 
— IV. This current saturation is not due to a series resistance 
between RGO and Ag, which is a common reason for forward cur- 
rent saturation in typical Schottky diodes, since the thick RGO diode 
with an identical ohmic Ag/RGO contact did not show a similar 
current saturation. Furthermore, the high bias up to 20 V in thick 
RGO did not show such an abrupt current saturation at all (see 
Supplementary Information, Section 2). In order to understand this 
peculiar phenomenon in the I-V characteristics of the thin RGO 
diode, an in-depth analysis of charge transport was performed on 
the basis of the I-V results. Considering the charge transport through 
an RGO/n-Si junction, tunneling (i.e., direct tunneling, DT, or 
Fowler- Nor dheim, FN) is highly possible for thin RGO/n-Si since 
the average thickness of RGO is on the order of several nanometers. 
The tunneling current and applied bias voltage are expressed by the 
following equations 21 : 
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Figure 1 | AFM images of (a) thick (—390 nm) and (b) thin (—5.7 nm) RGO layers, (c) Raman spectra for a thin RGO layer. 
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where m* is the effective mass of an electron, 0 is the barrier height, d 
is the thickness of RGO, e is the unit electrical charge, and H is 
Planck's constant. Typically, in tunneling current measurements, a 



transition from DT to FN tunneling takes place at a transition voltage 
(Vtrans) defined as V trans = &le. V trans can be found from either a plot 
of ln(l/V 2 ) versus ln(l/V) (i.e., a DT plot) or a plot of ln(l/V 2 ) versus 
(1/V) (i.e., FN plot). Figure 3b shows a forward current DT plot for 
thin and thick RGO/n-Si diodes. For thin RGO/n-Si diodes, an 
inflection point is observed, where the initial DT current (0 < V < 
0.4 V) transforms to an FN current at ~1 V. It has been established 
that in nanoscale junctions using carbon nanotubes 22 , 7i-conjugated 
molecules 23 , and Si nanosheets 21 , the charge transport undergoes 
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Figure 2 | XPS spectra for (a) C 1 s and (b) VB edge for thick and thin RGO and GO. (c) The O/C stoichiometric ratio and (d) the C 1 s bond fraction (C- 
C, C-O/C = O, and C-OOH) for thick and thin RGO and GO samples. 



consecutive transitions from (i) DT to (ii) FN to (iii) final breakdown. 
This trend is also observed for epitaxial graphene/graphene oxide 
junctions 20 . The DT plot shown in Fig. 3b follows the same charge 
transport behavior up to 1 V as reported in other studies; however, 
above 1 V, abrupt current saturation governed by second DT trans- 
port phase (i.e., DT (II)) is induced. Therefore, the thin RGO/n-Si 
diode has a striking multimode transition in charge transport — DT 

(I) — > FN — » DT (II) — disobeying the single charge transport trans- 
ition (DT — > FN) observed in typical nanoscale junction diodes. DT 

(II) occurs only in a thin RGO diode, implying a relation to the 
dimensional confinement that results in a dynamic tunneling barrier 
dependent on an electric field (E field). On the other hand, a thick 
RGO diode shows a significantly different transport mechanism that 
is more likely based on Schottky emission (SE) at different slopes at 



low and high bias respectively, as shown in the inset of Fig. 3b. The 
low linear slope of high bias (at >2 V) in SE plot might be due to the 
partial ohmic resistance. Other tunneling transport mechanisms are 
excluded for thick RGO because of its thickness of 390 nm. Based on 
the figure of merit for general tunneling mechanism, under assump- 
tion of same effective mass of electrons between two samples, tunnel- 
ing probability exponentially depends on thickness as expressed by 
exp(-t thick Rco/tthin rgo) equal to ~10~ 30 . Therefore, no tunneling 
in thick RGO is readily expected. This is confirmed in the inset of 
Fig. 3b and its transport mechanism is governed by Schottky emis- 
sion not by tunneling. 

A thin RGO diode shows more interesting changes in the I-V 
characteristics when a high bias voltage is applied. Figure 4 shows 
the I-V curve for a thin RGO/n-Si diode for bias voltages up to 10 V 
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Figure 3 | (a) A typical I-V curve on the logarithmic current scale and (b) DT plot for thin and thick RGO/p-Si diodes. The inset figures in (a) and (b) 

are the schematic illustrations of the RGO/n-Si diode structure and a Schottky- emission plot for a thick RGO diode, respectively. 
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under two different ambient conditions — air and N 2 . Under both 
ambient conditions, current saturation occurs at ~1 V. However, 
in air, the current rises again from the saturated current level above 
~4 V. In contrast, no current rise is observed in the N 2 environment. 
The DT plot indicated that charge transport in the current increase 
region (>4 V) is most likely due to DT transport. On the basis of the 
diode dependence on ambient conditions, two crucial conclusions 
are drawn: (i) the current saturation onset at 1 V is due to intrinsic 
properties of the diode, regardless of the atmosphere, and (ii) at a 
high bias voltage, electrically induced reactions between the surface 
of a thin RGO diode and gas species (only in air) can be considered as 
a strong modifier of the charge transport. The dependence of the 
conductivity change in a thin RGO diode on the ambient gas is 
significant: X 10 difference in the DT (II) region and X 100 difference 
in current rise region (>4 V). In thick RGO diode, same I-V curves 
were measured both in Air and N 2 ambient without any current 
saturation as in thin RGO (Data not shown here). 

A mechanism for the multimode transition of charge transport in 
the thin RGO/n-Si diode is proposed, as shown in Fig. 5a. We first 
consider the possible physical origin of the current saturation onset 
and the DT (II) region. From Fig. 3b, it is inferred that a bias-induced 
tunneling barrier is formed at onset voltage and this barrier is dis- 
tinguished from the tunneling barrier of the DT (I) region. Because 
this is related to the intrinsic properties of the diode, it is highly 
possible that negative O ion migration in RGO under an E field forms 
the new tunneling barrier for the DT (II) region. Recently, S. Kwon 
et al reported that metastable O ions in ultrathin Cu/single layer 
graphene junctions (on the scale of several nanometers) can be 
reversibly drifted by E field control, and also suggested on the basis 
of density functional calculations that reversible O migration can be 
accessed via metastable O sites on the graphene surface 16 . 
Interestingly, the reversible change in resistance due to O ion drift 
has also been observed at a bias voltage of ~ 1 V 15 , which is equal to 
the current saturation voltage in this study. Another theoretical 
study 24 predicted that gate voltage control of O diffusion is possible 
by the creation of metastable weak bonding from additional electrons 
supplied by the gate bias, thereby lowering the diffusion barrier 
energy (0.15 eV) by an order of magnitude from the measured value 
(1.5 eV) on the graphene surface in thermal equilibrium 25 . 
Therefore, due to the electrostatic interaction between E-field and 
O ions in RGO, the positive bias applied to the Au electrode on RGO 
surface will move negative O ions electrostatically toward the RGO/ 
Ag junction. This is possible because RGO has a significant amount 
of unsaturated metastable graphene sites available for O ion migra- 
tion. However, O ions distributed in thick RGO should migrate over 
a much longer distance by a factor of 68 than that in thin RGO in 
order to be reversibly localized at either the surface of RGO or RGO/ 



30 
25 
— 20 

1 

r 15 



10 



— I — I — I — 1 — 


r-r> — i — 1 — i — 1 — i — 1 — i — 1 — : 

A 

1 




1 

1 ' 


x100 A 


\ \ c 
» * 

- * * 


•fa 0 ...,: 


y »• 




In Air 

7 >nN o 

J 2 




1 1 1 1 


I i 1 i 1 i 


i i i i 



10- 



10 



* ^ 



10° 



10 



o 

9 o 



-4 -2 0 2 4 6 8 10 
Bias (V) 

Figure 4 | Current and conductivity vs. voltage plot for a thin RGO/p-Si 
diode in air and N 2 atmospheres. 



Si interface. This requires much higher activation energies of O ion 
migration due to the higher number of collisions and interactions 
with lattice atoms than those required in thin RGO. Therefore, this 
limits the effect of O ions for charge transport change or no DT II 
region occurs in I-V curve of thick RGO diode. 

As a result of the O ion drift, an O-rich region is formed at RGO/ 
Ag interfaces, acting as the transient DT (II) barrier formed at a bias 
of > 1 V. This DT (II) barrier formation/disruption occurs reversibly 
in a way that the backward bias (from high to low voltage) recovers 
the original I-V curve, as shown in a hysteresis I-V analysis 
(Supplementary Information, Section 2). As another supporting 
evidence for O migration under E-field, the thermal dependence of 
the I-V characteristics was evaluated in order to measure the thermal 
activation energy of O-diffusion under E-field condition, as shown in 
Fig. 5b. The substrate temperature is varied from room temperature 
to 100°C in vacuum (at 5 mTorr range) in order to minimize any 
possible surface reactions from ambient air to focus on charge con- 
duction mechanism under E-field driven O diffusion from RGO. In 
Fig. 4b, the DT (II) region is eliminated as increasing temperature 
above 50°C. A continuous current increase is observed without the 
apparent DT (II) region while the DT (I) region and reverse current 
are almost identical to those observed at room temperature. 
Therefore, this behavior serves as evidence that O ions near the 
surface, which are involved in the transient DT (II) barrier formation 
at room temperature, are thermally more activated under the high E- 
field condition. At temperatures higher than 50°C, thermal activa- 
tion of O ion diffusion and relevant carrier transport becomes more 
aggressive to result in the forward/reverse current increase at any 
voltage, thereby losing the rectifying junction characteristics. On the 
basis of the analysis of Fig. 5b, the temperature-dependent conduc- 
tance variation can be thought of as a thermal activation process 
related to OH and O diffusion and relevant carrier transport in 
RGO. Furthermore, the thermal conductance response is divided 
into two regions: (i) O-diffusion in low bias (i.e., DT (I) region) 
and (ii) O-diffusion in high bias (i.e., DT (II) region). Conductance 
values at 0.5 V and 2 V were considered according to the Arrhenius 
law as follows: 



Gccexp 



/c 5 T 



(2) 



where G is the conductivity, E A is the activation energy, k B is the 
Boltzmann constant, and T is the temperature. In this thermal I-V 
analysis, E A is related to thermally activated O diffusion and relevant 
carrier transport processes in RGO in the presence of an E field. 
Figure 5c shows the ln(G) versus l/k B T Arrhenius plot for the ther- 
mal dependence of the I-V curve and yields E A values of 0.17 and 
0.2 eV for the DT (I) and DT (II) regions, respectively. These activa- 
tion energy values are comparable to the metastable O diffusion 
barrier energy (0.15 eV) during electron supply, as predicted by 
calculations 24 . Although thermal reaction can result in complicated 
chemical/physical changes at such a Si/RGO nanoscale junction, 
only loosely bound O and OH in RGO can be thermally activated 
at the temperature range of 50-100°C under vacuum since all other 
bonds such as Si-Si, Si-O, and C-C have high bond strengths. 
Moreover, upon cooling of the diode measured at the elevated tem- 
perature in vacuum, I-V characteristic is clearly recovered to the 
original RGO property. This indicates that the thermal I-V change 
is not mainly due to thermal desorption of O from RGO, which is the 
irreversible chemical change. Indeed, O desorption in RGO has been 
also reported to occur aggressively at only high temperature 
>100°C 26 . We also should mention that thermal process at 
>100°C has already been gone through for RGO formation so its 
thermal stability can be assured. Then, the thermal change in I-V 
feature is mainly related to both (1) O migration in RGO and (2) 
thermally activated charge transport mechanism. Therefore, the 
activation energy from temperature-varied I-V curve should be also 
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related to (1) and (2), the major factor of which is not determined 
clearly at this point. However, O migration under E-field should be 
mainly or partly involved in the thermal activation process and the 
activation energy at 0.17 ~ 0.2 eV from Fig. 5b is comparable to the 
reported values of O diffusion barrier on graphene surface under E- 
field. Therefore, we regard that E-field can lower activation energy 
greatly in this study. Also, the wiggling I-V trace at the high voltage 
(>2 V) from >75°C indicates a typical transport character of 
reduced oxides with lots of defect sites (i.e., due to metastable O/ 
OH in RGO). Therefore, it is suggested that Fig. 5b supports that O 
ion loosely bound C is the main payer for charge transport mech- 
anism in RGO (OH is not a major species due to low bond faction as 
shown in Fig. 2d) and it is also controllable by external E-field due to 
the lowered activation energy extracted from Fig. 5b. 

Among possible mechanisms responsible for the second current 
rise in air following DT (II), RGO surface-gas reactions under an E 
field or during electron supply are considered the most probable 
cause. The constituent gas species in air are 0 2 , N 2 , and water. A 
reaction of water and N 2 is excluded on the basis of the results shown 
in Fig. 4. A highly possible electrochemical reaction of 0 2 and water 
at the RGO surface with a high density of electrons under a bias is 
expressed as follows: 

0 2 + 2H 2 0 + 4e~ -^40H " (3) 

Therefore, as long as this electrochemical reaction proceeds, OH ions 
can be supplied and adsorbed onto the RGO surface. The p-type 
electronic doping effect by the hydroxyl group on graphene has been 
widely studied 26,31 . If OH is chemisorbed on the RGO surface and acts 
as a p-doping enhancer, it will pull more electrons from carbon 
atoms, leading to a more positively charged RGO surface. Upon 
surface OH adsorption, the hole carrier density increases, and con- 
currently, the current flow by hole injection from RGO to n-Si 
rapidly increases. If such p-type doping is only limited at surface 



or near- surface regions by the OH diffusion length, the hole charge 
transport can be regarded as charge injection from the OH-rich 
surface to the RGO bulk until the collection of holes occurs at n-Si. 
This might lead to complex hole transport at the RGO/p-Si junction, 
as confirmed by the I-V analysis, revealing no convergence of the 
current at >4 V to a single charge transport mechanism. Therefore, 
it is suggested that the controlled multilevel conductivity of RGO is 
mainly based on a transient electrochemical effect, arising from an E- 
field- driven lattice O ion drift in bulk RGO and at the RGO/n-Si and 
RGO/metal interfaces. For electrochemical OH induced p-type dop- 
ing observed in thin RGO, this requires high flux of electrons with 
high kinetic energies to ensure over-potential to initiate electrochem- 
ical reaction. The number of high energy electrons (i.e., hot electrons) 
is much higher in thin RGO than in thick RGO due to tunneling 
process in the short distance less than carrier mean free path but in 
thick RGO, major portion of electrons is should be thermalized due 
to electron -lattice scattering (i.e., ion and phonon scattering). This 
explains why second current rise is not observed in thick RGO even 
upon applying high voltage up to —20 V. An RGO diode should be 
stable and reversible for practical applications such as chemical sen- 
sors and photodiodes. This aspect was studied on the basis of the 
hysteresis I-V analysis, showing a current deviation between forward 
and reverse sweeps only within the DT (II) region (Supplementary 
Information, Section 3). This is regarded as an indication of the 
electrochemically modulated current in the voltage region where 
OH formation/disruption according to Le Chatelier's principle and 
bidirectional lattice O ion drift occur in a forward/reverse sweep. 
Therefore, this is analogous to the statistically distributed O/metal 
ion drift- controlled resistive change in bistable switching oxides/ 
ionic materials 27 . In repeated hysteresis I-V measurements, biasing 
at 1 V acts as a "reset" voltage for recovery to an identical current 
level, regardless of the bias sweep direction., If the high- and low- 
current level in the DT (II) region becomes larger by some orders of 
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magnitude by device engineering, the hysteresis behavior provides a 
potential for use as resistive memory. The peculiar variation of the 
conductivity in this study is considered to be based on a transient 
electrochemical effect consisting of (i) an E-field-driven lattice O ion 
drift and (ii) surface OH p-type doping. It should be noted that a 
similar mechanism for electrochemical reduction/oxidation of GO 
was proposed such that each half of the reaction between GO and 
H 2 0 forms H + ions at the cathode and OH" ions at the anode, 
respectively 28 . These are very interesting and special properties of 
the RGO/Si diode geometry under the out-of-plane charge transport 
mode in nanometer-scale thin RGO and provide an important basis 
for the application of O bistability. This finding can be potentially 
utilized for RGO -based ionic devices such as chemical sensors in 
OH-related chemical environment (e.g., pH sensor) and 
electrocatalysis. 

Methods 

Expanded graphite was prepared via microwave-assisted thermal expansion of highly 
ordered pyrolytic graphite (SPI-2 grade from SPI). A graphene oxide dispersion was 
synthesized from the expanded graphite using the modified Hummers method and 
was reduced by the chemical reduction method. A dispersion of GO-hydrazine was 
prepared by mixing of as -synthesized GO with hydrazine monohydrate, followed by 
dilution (150 times) with a mixture of DI water and ethanol and sonication for 1 min. 
The as-prepared RGO dispersion was spray-coated onto a preheated (240°C) n-type 
doped and chemically cleaned Si wafer and onto a Si wafer with a thermally grown 
300-nm-thick oxide layer (Si/Si0 2 wafer), respectively. The resultant RGO thickness 
was in the range of ~5 nm for the thin layer and several hundred nanometers for the 
thick layer, respectively, as measured by step height analysis using atomic force 
microscopy. Raman spectra were recorded in the backscattering configuration using a 
micro-Raman spectrometer (Labram, ISA Groupe Horiba), equipped with a He-Ne 
laser (532 nm, —0.8 mW at the sample). For diode fabrication, the top Ag (100 nm) 
electrode on RGO was formed by e-beam evaporation and the rear electrode on the p- 
Si wafer was produced by applying a Ag paste. The chemical bonds of RGO were 
analyzed using XPS [ESCA Lab-22201, VG] with an Mg source. The binding energy of 
each element was calibrated using C-C bonds (284.5 eV) in the C 1 s binding state. 
The electrical characteristics of a diode and a resistor were measured by a vacuum 
chamber current-voltage (I-V) probe station connected to a Keithley 4200 semi- 
conductor analyzer. I-V analyses were carried out under diverse atmospheres such as 
vacuum (10~ 3 torr), air, and N 2 in order to investigate electrochemical changes of the 
RGO surface. 
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